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We analyze the azimuthal dependence of the heavy-quark-initiated contributions to the lepton- 
nucleon deep inelastic scattering (DIS). First we derive the relations between the parton level semi- 
inclusive structure functions and the helicity 7*Q cross sections in the case of arbitrary values of the 
heavy quark mass. Then the azimuth-dependent 0{as) lepton-quark DIS is calculated in the helicity 
basis. Finally, we investigate numerically the properties of the cos ip and cos 2ifi distributions caused 
by the photon-quark scattering (QS) contribution. It turns out that, contrary to the basic photon- 
gluon fusion (GF) component, the QS mechanism is practically cos 2(y3- independent. This fact implies 
that measurements of the azimuthal distributions in charm leptoproduction could directly probe the 
charm density in the proton. 
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I. INTRODUCTION 



The notion of the intrinsic charm (IC) content of the proton has been introduced over 25 years ago in Refs P, 0]. 
It was shown that, in the light-cone Fock space picture [3,13], it is natural to expect a five-quark state contribution, 
] \uudcc), to the proton wave function. This component can be generated by gg cc fluctuations inside the proton 
where the gluons are coupled to different valence quarks. The original concept of the charm density in the proton 
1,0] has nonperturbative nature since a five-quark contribution \uudcc) scales as where m is the c-quark mass 



> 

^\ ' A decade ago another point of view on the charm content of the proton has been proposed in the framework of the 
, variable flavor number scheme (VFNS) @, III- Within the VFNS, the mass logarithms of the type In (Q^/m^) are 
resummed through the all orders into a heavy quark density which evolves with according to the standard DGLAP 
, evolution equation. Hence this approach introduces the parton distribution functions (PDFs) for the heavy quarks 
■ and changes the number of active flavors by one unit when a heavy quark threshold is crossed. Note also that the 
C^i charm density arises within the VFNS perturbatively via the g —> cc evolution. Some recent developments concerning 
the VFNS are presented in Refs. [1, i, [13, [O, [H ■ 
] Presently, both nonperturbative IC and perturbative charm density are widely used for a phenomenological descrip- 
t tion of available data. (A recent review of the theory and experimental constraints on the charm quark distribution 
' can be found in Refs. [Ij, [l?]). In particular, practically all the recent versions of the CTEQ [l^ and MRST [l6| sets 
• ^ , of PDFs are based on the VFN schemes and contain a charm density. At the same time, the key question remains 
' open: How to measure the charm content of the proton? The basic theoretical problem is that radiative corrections 
^ to the leading order (LO) predictions for the heavy quark production cross sections are large: they increase the Born 
- - 1 level results by approximately a factor of two at energies of the fixed target experiments. On the other hand, per- 
turbative instability leads to a high sensitivity of the theoretical calculations to standard uncertainties in the input 
QCD parameters: to, /i^, iip, Aqcd and PDFs. For this reason, one can only estimate the order of magnitude of the 
pQCD predictions for the heavy flavor production cross sections [17„ ,1^. 

At not ver y high e nergies, the main reason for large NLO cross sections of heavy flavor production in ^g [H,[13], 7*ff 
and gg |22l.T23l. l23. l25l] collisions is the so-called threshold (or soft-gluon) enhancement. This strong logarithmic 
enhancement has universal nature in the perturbation theory since it originates from incomplete cancellation of the 
soft and coUinear singularities between the loop and the bremsstrahlung contributions. Large leading and next-to- 
leading threshold logarithms can be resummed to all orders of perturbative expansion using the appropriate evolution 
equations [1^, [13, [13] ■ Soft gluon resummation of the threshold Sudakov logarithms indicates that the higher-order 
contributions to the heavy flavor production are also sizeable. (For a review see Refs. (29l. [sol. [3l| ) . 
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Since production cross sections are not perturbatively stable, it is of special interest to study those observables that 
are well-defined in pQCD. A nontrivial example of such an observable was proposed in Refs. [H, [H, [13, [sH] where the 
azimuthal cos2(^ asymmetry in heavy quark photo- and leptoproduction has been analyzed ^. In particular, the Born 
level results have been considered 132. 13411 and the NLO soft-gluon corrections to the basic mechanism, photon-gluon 
fusion (GF), have been calculated [33, [3J]. It was shown that, contrary to the production cross sections, the cos 2^3 
asymmetry in heavy flavor photo- and leptoproduction is quantitatively well defined in pQCD: the contribution of 
the dominant GF mechanism to the asymmetry is stable, both parametrically and perturbatively. This fact provides 
the motivation for investigation of the photon- (heavy) quark scattering (QS) contribution to the (ys-dependent lepton- 
hadron deep inelastic scattering (DIS). 

In the present paper, we calculate the azimuthal dependence of the next-to-leading order (NLO) ©(aem^s) heavy- 
fiavor-initiated contributions to DIS. To our knowledge, pQCD predictions for the (^-dependent cross sections in 
the case of arbitrary values of the heavy quark mass m and are not available in the literature. Moreover, there is 
a confusion among the existing results for azimuth-independent cross sections. 

The NLO corrections to the (/j-independent lepton-quark DIS have been calculated (for the first time) a long time 
ago in Ref. [1^ , and have been re-calculated recently in [13 • The authors of Ref. [13| conclude that there are errors 
in the NLO expression for a^'^^ given in Ref. [1^ ^ . We disagree with this conclusion. It will be shown below that a 
correct interpretation of the notations for the production cross sections used in [36} leads to a complete agreement 
between the results presented in Refs. [11], [1^ and present paper. 

As to the (^-dependent cross sections, our main result can be formulated as follows. Contrary to the basic GF 
component, the QS mechanism is practically cos 2(/3- independent. This is due to the fact that the QS contribution 
to the cos 2(/3 asymmetry is absent (for the kinematic reason) at LO and is negligibly small (of the order of 1%) at 
NLO. This fact indicates that the azimuthal distributions in charm leptoproduction could be a good probe of the 
charm density in the proton. In detail, the possibility of measuring the charm content of the proton using the cos 2lp 
asymmetry will be investigated in a forthcoming publication [soj . 

Concerning the experimental aspects, azimuthal asymmetries in charm leptoproduction can, in principle, be mea- 
sured in the COMPASS experiment at CERN, as well as in future studies at the proposed cRHIC [40, |4l| and LHeC 
[i^ l colliders at BNL and CERN, correspondingly. 

The outline of this paper is as follows. In Section |lTl we derive the relations between the parton level semi- 
inclusive structure functions and the helicity cross sections in the case of arbitrary values of the heavy quark 
mass. As explained in Ref. [i^, in the presence of non-zero masses, it is the helicity basis that provides the simplest 
connections between the hadron- and parton-level production cross sections. In Section IIIIl we present the NLO 
0{aemCis) predictions for the (/3-dependent lepton-quark DIS in the helicity basis. Our calculations are compared with 
available results in Section [Tvl In Section |Vl a numerical investigation of the costp and cos 2(^9 distributions caused by 
the QS contribution is given. In particular, we provide a simple parton level proof of the fact that the QS mechanism 
is practically cos 2(p- independent. Our conclusions are presented in Section [VII 



II. AZIMUTH-DEPENDENT STRUCTURE FUNCTIONS IN THE HELICITY BASIS 

In this Section, the helicity formalism for the semi-inclusive cross sections in the case of arbitrary values of 
the heavy quark mass is presented. This is a purely kinematical analysis, which will set the notation to be used later 
on. In fact, we extend the helicity approach proposed in Ref. [43] to the case of (^-dependent leptoproduction using 
the method formulated in Ref. [44} . 

We consider the semi-inclusive deep inelastic lepton-quark scattering. The momentum assignment will be denoted 

as 

m + Q{k) ^ l\i -q) + Q'{p) + X{px). (1) 
The following definition of partonic kinematic variables is used: 

^^FI' ^^2^' ^=Q^' (2) 



^ The well-known examples are the shapes of differential cross sections of heavy flavor production which are sufficiently stable under 

radiative corrections. 
^ For more details see PhD thesis [3^ . pp. 158-160. 
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FIG. 1; Definition of the azimuthal angle ip in the target rest frame. 



The differential cross section of tfie reaction ([T]), d^aiq, is defined in terms of the quark tensor Wq' 



^0 ,3 



dydQ^dip ^ y dzdQ^dip ~ {e ■ k)Q^ '"^ ^ (27r)32po '' 



(3) 



where = — is the 4-niomentum of the final lepton. In the target rest frame, the azimuth ip is the angle 
between the lepton scattering plane and the heavy quark production plane, defined by the exchanged photon and the 
detected quark Q' (see Fig. [T]). The covariant definition of ip is 



cos (/3 



QVl + 4Az2 
sm (p — = n ■ I, 



The explicit expression for the lepton tensor is: 



(4) 



(5) 



where ^ denotes a sum over all final helicity states and an averaging over all initial spin variables. The semi- inclusive 

spin 

quark tensor Wq" is defined as follows: 



W^-'{q,k,p)^l- {k\J^\p,Px){2nfd(^^Hq + k-p-px){px,p\J^''\k) 



Aw 



(6) 



X(px),spm 



where sums and integrals over all the unobserved final states X of momentum p'^ are implied. 

To construct the parton tensor describing the semi-inclusive 7**5 DIS, it is convenient to introduce two 4-vectors: 



q ■ k 



(7) 



{q- k)V^VTT4\^' ' Q 

that obey the following conditions: v ■ k ^ v ■ q ^ 0, v'^ = —1 and w-q = 0,w-k = w^— m? (l -I- 4Az2) /AXz^. In 
terms of w'^ and w'^ , has the following structure: 



q'^'^qv 



m 



(8) 



that obeys all the necessary conservation laws. In particular, Wq^q^ — 0. The scalar coefficients Wi (i = 1,2,/, ^4) 
are the semi-inclusive parton-level structure functions for the process ((TJ, Wi = Wi (^z,X,p^,q -p). 
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For parton-model considerations, is it convenient to use the so-called coUiner frames where the S-momenta of the 
virtual photon and initial quark are antiparallel to each other, q \\ (— fc). Evidently, an arbitrary coUiner frame can 
be obtained from the initial quark rest system with the help of a Lorcntz boost along q. Pointing the z-axis along g, 
we will have in a coUiner frame: 

= (0,«j_,0), ^i = rn' (9) 



gVl + 4Az2 ^ 1 



e^', e;; = -(|gl,Oi,go), (10) 



where Q = \J —q^ and Cq describes the longitudinal polarization of the virtual photon, 7*. It is also useful to define 
the scalar, e^, and transverse, e^, polarization vectors: 

e^-^, e^ = -^(0,Tl,-i,0). (11) 



Note the completeness relation 



e^l^e- + ete"* - e^e^* + e^e^* = (12) 



and the normalization for the physical states: 

e, -e: = (-1)^<5,„ (r,s = 0,±l). (13) 

One can see from Eqs. (I9ll0p that it is merely the scalar coefficient functions Wi {i = 1,2, 1, A) depend on the final 
quark momentum p in a coUinear frame. For this reason, we can integrate the semi-inclusive quark tensor Wq" {q, k,p) 
over p and obtain the inclisive quantity Wq^^q, k): 



The inclusive coefficient functions Wi{z, A) (i = 1, 2, /, A) are related to the semi-inclusive ones as follows: 

Wi{z,X) = I (^1^ (^1 + ^^^) Wi{z,X) = I j^^^Wj{z,X,p%,q-p),{15) 

Integrating WQ'{q,k) over the lepton azimuth ip defined by Eqs. one can reproduce the well-known expression 
for totally inclusive DIS: 

^ l'^,yW^'^(q,k) = -(^g^'^-^^W,iz,Q') + ^W2iz,Q'). (16) 

Note that the above relation can easily be obtained from Eq. (1141) taking into account that J^^d(pv'^v'^ — 

Now the cross section for the inclusive azimuth-dependent lepton-quark DIS can be written as 

d^o-ZQ _ y d^CTiQ _ y d^(T/Q _ al^y 



dzdQ^dip ~ z dydQ^d ip~ dzdyd(p ~ ^ 



(17) 



To derive the relations between the invariant and helicity structure functions, we use the completeness (|12p which 
implies that 
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where the quark and lepton heUcity structure functions (Frs and p^s, respectively) are defined as 

Frs - eJfW^Q.^.er, Prs = (-l)''+^erip.er. (19) 

Choosing the x-axis along v± defined by Eq. we obtain for the quark helicity structure functions Frs{z, A): 



F++ = F. 



F+o — Pq- 



F,o = -m + ^^T^w,, 



4Az2 



The lepton tensor prs {i, f) has the following form in the helicity basis; 



P++ = p — = 



P+O = Po + 



1 



-P- 



-Po- 



1-e 



Poo 
P+- 



2Q^i 
-P-+ -- 



(20) 



(21) 



l-e 



The quantity e measures the degree of the longitudinal polarization of the virtual photon in the Breit frame [4J] . The 
covariant definition is: 



2(l-y- Az2y2 



(22) 



l + (l-y)2 + 2Az2y2- 

In terms of the helicity structure functions, the azimuth-dependent inclusive Icpton-quark cross section has the form 



where 



dzdQMip 1-i 



Ft — F++, 



Friz, X) + sFLiz, X) + eFAiz, X) cos 2ip + 2 v/e(l + s) Fi{z, A) cos (p 



Fa = -F+^^ 



Fi = F^o- 



(23) 



(24) 



Likewise, using Eqs. (|18ll22p . one can easily express the semi-inclusive IQ cross section defined by Eq. ^ in terms of 
the corresponding helicity structure functions Frs (z, X, k ■ p, q ■ p) = e!^{q, k)WQ,fj_i,{q, k,p)e'^*{q, k). 
Sometimes, instead of the structure functions Frs, the helicity 7*Q cross sections are used: 



(7i(z, A) = 



Sn'^aem z 
Q2V1 + 4Az^ 



F,(z,A), {i^T,L,A,I), 



(25) 



where ax — <5'++, o'l = <Tooj = and aj — (t_q. Since y <C 1 in most of the experimentally reachable 

kinematic range, it is the the quantities (T2 and F2 that can effectively be measured in (^-independent DIS: 



0-2(2, A) = aT{z,X) + o-l(z,A), 



F2izA) 



2z 



1 + 4Az2 

In terms of the quantities (jj , the cross section of the reaction ([T]) can be written as 
aem Vl + 4Az2 



FTiz,X)+FLiz,X) 



(26) 



d^ 



<^IQ 

dzdQ^dtp ~ (27r)2 zQ^ 1 - e 



0'2 



(z. A) — (1 — £)(7l(z, A) + e<TAiz, A) cos2ip + 2y^e{l + e) &i{z, X) cos tp . (27) 



In Eqs. (|26|27p . ctt (ctl) is the usual 7* cross section describing heavy quark production by a transverse (longitudinal) 
virtual photon. The third cross section, a a, comes about from interference between transverse states and is responsible 
for the cos 2ip asymmetry which occurs in real photoproduction using linearly polarized photons [s^. [ssl . [3a | . The 
fourth cross section, tr/, originates from interference between longitudinal and transverse components [44| . 
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FIG. 2: The LO (a) and NLO (b and c) photon-quark scattering diagrams. 



III. PHOTON-QUARK SCATTERING CROSS SECTIONS AT NLO 

At leading order, 0{aem), the only quark scattering subprocess is 

l*iq) + QikQ)~^QiPQ). (28) 
The 7*Q cross sections, af^"^ (i — 2,L,A,I), corresponding to the Born diagram (see Fig. [2^) are: 



a^°^(z,A) = (7i3(z)v/l + 4Az2^(l-^), 



(z,A) = dsiz) \ , ^(l-z), (29) 
V 1 + 4Az^ 

ai°)(^,A) - (z,A) = 0, 



with 



em 



where eg is the quark charge in units of electromagnetic coupling constant. 

To take into account the NLO 0{aemCts) contributions, one needs to calculate the virtual corrections to the Born 
process (given in Fig.[2t) as well as the real gluon emission (see Fig.[2l3): 

l*{q) + Q{kQ)^Q{PQ) + g{Pg). (31) 

The NLO (/j-dependent cross sections, (t^^^ and (t\^\ are described by the real gluon emission only. Corresponding 

contributions are free of any type of singularities and the quantities ct^^'' and (T^^^' can be calculated directly in four 
dimensions. 

In the i^-independent case, (jg^'' and ^^]^\ we also work in four dimensions. The virtual contribution (Fig. 
contains ultraviolet (UV) singularity that is removed using the on-mass-shell regularization scheme. In particular, we 
calculate the absorptive part of the Feynman diagram which has no UV divergences. The real part is then obtained 
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by using the appropriate dispersion relations. As to the infrared (IR) singularity, it is regularized with the help of an 
infinitesimal gluon mass. This IR divergence is cancelled when we add the bremsstrahlung contribution (Fig. [SJd). 
The final (real+ virtual) results for cross sections can be cast into the following form: 



A) = ^Ci.CTi3(l)Vl + 4A(5(l - z)f-2 + 41nA - VTTAX In 



1 + 2A 



21nD(z,A) 



Vl + 4Az2 



2Li2(r2)+4Li2(-r) 

Vl + 4A L 

+31n^(r) -41nr + 41nrln(l + 4A) -21nrlnA | 
^1 - 3z - 4z2 + Qz^ + Sz" - 
+6Az (3 - 18z + 13^2 + lOz^ - 8z'') 
+A\^z^ (8 - 77z + 65z2 - 22^) 

+ 16A^z3 (1 - 21z + 12z2) - 128A''z5 

(1 + z + 2z2 + 2z^) + 2Az (2 - llz - llz^) + SA^z^ (1 - 9z) 
8(1 + 4A)2z4 4(1 + 2A)(1 + AXfz^ \xiD{z, A) 



d'i\z,X)^^CF<TB{l) 



2A 



VI + 4A 



^(l-z){- 



-2 + 4 In A 



4A 



Inr 



VI + 4Az2 
1 + 2A 



(1-^)h 



-^^'^^^'\l + L^f/2 I [1 - (1 - A)z]^ 



, , 2Li2(r^) + 4Li2(-r) 

ViTlA VTTIaL ^ ' 

•i\n{r) - 41nr + 41nrln(l + 4A) - 21nrlnA | 
(1 - zf ~ Az (13 - 19z - 2z2 + 8z3) 

-2\^z^ (31 - 39z + 8z2) 
-8A='z=' (10 - 7z) - 32AV^ 



2\z'^\iiD{z,\) 



[3 + 3z + 16Az] 



VI + 4Az2 

8A(1 + 4A)z4 4A(1 + 2A)(1 + 4A)z4 lni:i(z, A) 



(1-^)+ 



as 



a^^\z,X) = ^CpcfBiz) 



z{i-z) r 



VI + 4Az2 

2AlnL>(z, A) 



27r 



a«(z,A) 



8V2 



Cfctb{z] 



(l + 4Az2)2 [i_(i_A)zf/' 



-(1 - z)(l + 2z) - 4Az (10 - lOz - z^ + 2z^) 



(32) 



(33) 



-8A^ 



25 - 29z + 8z2) - 96A^z3 (3 - 2z) - 128A''z'' (35) 



+8v/Az [1 - (1 - A)z] [1 - + Az(13 - llz) + A\^z^{7 - 4z) + IGA^z^ 



In Eqs. ([5355]) . Cf = (A^c " l)/(2A^c), where is number of colors, while 



D{z,\) = 



1 + 2Az - VT+4Ai2 



1 + 2Az + Vl + 4Az2 ' 
The so-called "plus" distributions arc defined by 



= V7^(z = l,A) = 



V1 + 4A- 1 



Vl + 4A + r 



[5(z)]+=ff(z)-5(l-z) 1 dC5(C). 



(36) 



(37) 
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For any sufficiently regular test function h{z), Eq. ((57)) gives 



dz h{z) 



'ln''(l-z) 
1 - z 



, ln'^(l — , ,^ , ln'^"'""^(l — a) 
dz — ^ [h{z) - h{l)] + h{l) ^ ^ 



1 



1 



(38) 



IV. COMPARISON WITH AVAILABLE RESULTS 



For the first time, the NLO 0{aemCts) corrections to the iy9-independent IC contribution have been calculated a long 
time ago by Hoffmann and Moore (HM) [36]. However, authors of Ref. [36*1 don't give explicitly their definition of the 
partonic cross sections that leads to a confusion in interpretation of the original HM results. To clarify the situation, 
we need first to derive the relation between the lepton-quark DIS cross section, daig, and the partonic cross sections, 
cr*^^-' and cr'-^'', used in [1^. Using Eqs. (C.l) and (C.5) in Ref. ^36], one can express the HM tensor cr'^ in terms of 
"our" cross sections (T2 and ctl defined by Eq. ([77]) in the present paper. Comparing the obtained results with the 
corresponding definition of a'^ via the HM cross sections ct*^^^ and cr'^^ (given by Eqs. (C.16) and (C.17) in Ref. 36]), 
we find that 



a2{z,X) = <7s(z)Vl + 4Az2(7(2)(z,A), (39) 

a.izA) ^ -== p^'(z,A) + 2AzV2)(z,A)l . (40) 
Vl + 4Az2 L J 

Now we are able to compare our results with original HM ones. It is easy to see that the LO cross sections (defined 
by Eqs. (37) in and Eqs. (^5)) in our paper) obey both above identities. Comparing with each other the quantities 

a^^ and a2^^ (given by Eq. (51) in [sl] and Eq. ([5^ in this paper, respectively), we find that identity ([5^ is satisfied 
at NLO too. The situation with longitudinal cross sections is more complicated. We have uncovered two misprints in 
the NLO expression for a^^'^ given by Eq. (52) in First, the r.h.s. of this Eg. must be multiplied by z. Second, 
the sign in front of the last term (proportional to S{1 — z)) in Eq. (52) in Ref. [36( must be changed ^. Taking into 
account these typos, we find that relation (PD|) holds at NLO as well. So, our calculations of (72 and ctl agree with 
the HM results. 

Recently, the heavy quark initiated contributions to the iy9-independent DIS structure functions, F2 and F^^ have 
been calculated by Kretzer and Schienbein (KS) [37]. The final KS results are expressed in terms of the parton level 
structure functions Hf and Using the definition of Hf and given by Eqs. (7,8) in Ref. [s^l, we obtain that 



"^^'^ = ^ Vl + 4A Vl + 4Az2' = ^""^^VtTIa;^ ' ^''^ 



where ctt = o'2 — itl and ctl are defined by Eq. ((?7|) in our paper and ^' — z (l + Vl + 4A) /(l + vT+4Az^) . To 
test identities (PT|) . one needs only to rewrite the NLO expressions for the functions Hfi^', A) and -ff|(^'. A) (given in 
Appendix C in Ref. [37*1) in terms of variables z and A. Our analysis shows that relations (|jT|l hold at both LO and 
NLO. Hence we coincide with the KS predictions for the 7*Q cross sections. 

However, we disagree with the conclusion of Refs. [13, [la] that there are errors in the NLO expression for cr'^^^ 
given in Ref. {s^. As explained above, a correct interpretation of the quantities cr'^^^ and cr^^' used in [3^] leads to a 
complete agreement between the HM, KS and our results for (^-independent cross sections. 

As to the ((f-dependent DIS, pQCD predictions for the j*Q cross sections aA{z, A) and (7/(z, A) in the case of arbitrary 
values of and are not, to our knowledge, available in the literature. For this reason, we have performed several 
cross checks of our results against well known calculations in two limits: and 0. In particular, in 

the chiral limit, we reproduce the original results of Georgi and Politzer [i^ and Mendez [i^ for ai{z, A 0) and 
(Ta{z,X 0). In the case of 0, our predictions for a2{s,Q^ 0) and (74(5,0^ 0) given by Eqs. p2|34p 



reduce to the QED textbook results for the Compton scattering of polarized photons [47 1 



^ Note that this term originates from virtual corrections and the virtual part of the longitudinal cross section given by Eq. (39) in Ref. [3^ 
also has wrong sign. 
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FIG. 3: c^'"' (r;, A) and cf'"\ri,X) coefficient functions at several values of A. 
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V. SOME PROPERTIES OF THE AZIMUTH-DEPENDENT CROSS SECTIONS 



To perform a numerical investigation of the inclusive partonic cross sections, ai (i — 2, L, A, I), it is convenient to 



introduce the dimensionless coefficient functions c. 



(n,0 



z — / z — / \ jjj^ 

n=0 /=0 ^ 



(42) 



where /j, is a factorization scale (we use fi = i^ip = Hr) and the variable rj measures the distance to the partonic 
threshold: 



s 1 — z 

'n^ — -'^ = — > 



{q + kf 



(43) 



Our analysis of the quantity c'^'^\ri, A) is given in Fig. [21 One can see that c^'^"* is negative at low (A^^ < 1) 
and positive at high (A~^ > 20). For the intermediate values of Q^, c^2'°\ri, A) is an alternating function of rj. 
Let us discuss the coefficient function c^'°^(s, Q^) for the case of on-mass-shell photon, 0. In this hmit. 



MO) 



is, 



0) = SttCj 



(s — m^)2 



s + 7V? , m? "I 
2 H ^ In — 



(44) 



Considering now the threshold behavior of Eq. ((44|) . we find: lim c)^' {s,Q ^ 0) = — 47rCF/3. Taking also into 

s — *ni^ 

account that lim c"^'^\s^Q'^ ^ 0) = 0, we see that the mass-shell, — > 0, and threshold, s —t m?, limits do not 

s — 

commutate with each other for the quantity c'^''^\s , Q"^) . This property of the cross section c'^''^\s , Q"^) illustrates 
the well-known fact that there is no, generally speaking, a smooth transition between the lepto- and photoproduction. 

Our results for the coefficient function Cj°'°^ (77, A) at several values of A are presented in Fig. [3] It is seen that c^°'"' 
is negative at all values of rj and A. Note also the threshold behavior of the coefficient function: 



4A 



0(77). 



(45) 



This quantity takes its minimum value at Am = 1/4: cf'^\r] = 0, A„i) = —tt'^Cf/ (2^2)- 
In the chiral limit, rri^ 0, the (/9-dependent cross section are as follows: 



■'f-°\z, A ^ 0) = -^n'CFXzil + 2z)J—— + 0[\'). (46) 
4 y \ ~ z 



V2 



Let us analyze the numerical significance of the cos 95- and cos 2(^-distributions for the QS component. It is difficult 
to compare directly the a''^{z,X) and a'f'\z,X) cross section given by the usual functions and ((55)) with the 
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FIG. 4: The quantities a'x\N , X) / a!^''' {N , \) {left panel) and 2^/2 (iV, A)/o-^°' (iV, A) {right panel) at several values of A. 



99-independent contributions a^\z,X) and a^\z,X) described by the generahzed functions ([2^ and (j32|) . For this 
reason, we consider the Melhn moments of the corresponding quantities defined as 

1 

a,(7V,A) =y"a,(z,A)z^-Mz, {i^2,L,A,I). (47) 



The Melhn transform of the Born level cross sections is trivial: a^\N,X) = o's(l) Vl + 4A. The Mellin moments of 
the NLO results have been calculated numerically. We use for as{i^F) the one-loop approximation with A4 = 326 
MeV, = + and m = 1.3 GeV. 

The left panel of Fig. [4] presents the ratio a'^^\N , X) / u'"2\^ ^) as a function of N for several values of variable 
A: A^^ — 1,4,10,20 and 100. One can see that this ratio is negligibly small (of the order of 1%). Moreover, our 

analysis shows that the ratio a'"^ [N, X)/a^\N, A) is less than 1.5% for all values of A and > 0. This implies that 
the photon-quark scattering contribution is practically cos 2(/3-independent. 

In the right panel of Fig. [H the iV-dependence of the ratio 2y/2 a''p {N , X) / a^2'\N , X) is given for the same values 
of A. One can see that this ratio is of the order of 10-15% at small N and sufficiently high Q^. This fact indicates 
that the cos (/^-distribution caused by the QS component may be sizable. 



VI. CONCLUSION 



We conclude by summarizing our main observations. In the present paper, we have studied the azimuth-dependent 
photon- (heavy) quark DIS at NLO. It turns out that the cos2(^ dependence of the QS mechanism is negligible while 
the cos ip one may be sizable. The situation is diametrically opposite to the one that takes place for the basic GF 
contribution. It is well known that the GF predictions for the azimuthal cos2(/? asymmetry in heavy quark photo- 
[Hjlill and leptoproduction [H, [H, [HO] are large (about 20%). As to the costp dependence of the GF contribution, 
it vanishes at LO due to the charge symmetry Q Q [5l| . 

Since the GF and QS mechanisms have strongly different azimuthal distributions, one could expect that measure- 
ments of the (^-dependent DIS will directly probe the charm content of the proton. In detail, hadron level predictions 
for the azimuthal asymmetries as well as the possibility to discriminate experimentally between the GF and QS 
contributions will be investigated in Ref. [soj . 
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